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Abstract 

Experiments are carried out with a globally coupled, externally forced population of limit-cycle 
electrochemical oscillators with an approximately unimodal distribution of heterogeneities. Global 
coupling induces mutually entrained (at frequency ui) states; periodic forcing produces forced- 
entrained (wp) states. As a result of the interaction of mutual and forced entrainment, resonant 
cluster states occur with equal spacing of frequencies that have discretized frequencies following a 
resonance rule u;„ = nui — (n — 1)u;f- Resonance clustering requires an optimal, intermediate global 
coupling strength; at weak coupling the clusters have smaller sizes and do not strictly follow the 
resonance rule, while at strong coupling the population behaves similar to a single, giant oscillator. 
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I. INTRODUCTION 

Rhythms, often generated as synchronization of oscillator populations can also be 
exposed to external global forcing or feedback. The resulting dynamical behavior depends on 
the types of oscillators (smooth, relaxation, chaotic), topology of interactions (local, global, 
or network), and on the magnitude of forcing frequencies relative to the inherent frequencies 
of the oscillators. A plethora of interesting behavior has been observed with forcing in 
reaction-diffusion systems or locally coupled, identical oscillators: some examples include 
labyrinthine standing waves [2|, standing wave patternsis], resonant phase patterns 4\, Bloch- 
front turbulence js I, localized clustersjsl, spirals with hypocycloidal trajectoriesj?], wave traps 



and twisted spirals [8|, and pacemaker entrainments[9||. 

Less attention has, however, been paid to the analysis of heterogeneous populations of 
oscillators with a distribution of natural frequencies. A simple example is a globally coupled 



population with global forcing. Following Kuramoto's heuristic argument [10|, it is expected 
that global coupling would produce a mutually entrained state; this state behaves as a giant 
oscillator and becomes entrained to the frequency of an external forcing signal. However, 
■;he frequency adjustment process of the mutually entrained state is not trivial: Sakaguchi 



Ufl studied the effects of the external fields on mutual entrainment by analysis of a phase 



model with unimodal heterogeneities. The simulation results show that a transition from 
the mutual entrainment to the forced entrainment occurs as forcing strength is enhanced; in 
between, two and more plateaus are seen in the frequency-vs-natural-frequency plots which 
indicate the formation of multiple, resonance clusters. Because of the interaction of the 
mutually entrained cluster (with frequency of Ui) and the forced cluster (cuf), new resonant 
clusters with frequency of u;„ = nui — {n— l)ujf are formed at weak forcing strengths. Note 
that these discretized frequencies are equally spaced with /S.u = Ui — Uf. It was also shown 
that the order parameter exhibits large amplitude oscillations when two major clusters are 



formed. 
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In this paper we investigate experimentally resonance clustering in a chemical system 
(the oscillatory electrodissolution of a nickel electrode array in sulfuric acid solution) with 
global coupling and forcing with limit cycle oscillators. The mutual and forced entrainment 
states are identified, and their interactions are analyzed. Features of resonance clustering 
are compared as coupling strength and forcing amplitude and frequency are varied. Numer- 
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ical studies are carried out to confirm the experimental finding of resonance clustering in 
ordinary difi^erential equation models and to investigate the features of the dynamics in a 
large parameter space. 



II. EXPERIMENTS 

A. Experimental Setup 

The experiments were carried out in a standard three electrode electrochemical cell con- 
taining 3 mol/dm^ sulfuric acid at 11 °C with Ni working, a Hg/Hg2S04/K2S04 reference, 
and a Pt counter electrodes [3, The currents {ik{t)) of Ni electrodes (64 1-mm diameter 
electrodes in an 8 x 8 geometry with 2 mm spacing) were measured at 100 Hz. The potential 
of each electrode was held at potential V versus Hg/Hg2S04/concentratedK2S04 reference 
electrode. The external forcing is added to the applied potential V{t) = Vq + fein(27r/Ft). 
Vo = 1.09 V. The electrodes were connected to the potentiostat through one series (collec- 
tive) resistor, Rg, and through 64 parallel resistors jlj, Q. The interaction strength K was 
controlled through the external resistors K = (i?s/i?tot)/(l — Rs/Rtot) where i?tot = lO.lil 
was the total resistance. 

The phases and frequencies of the oscillators were determined with the Hilbert transform 
method 15 1 from time series data of current ik{t) jl3|. An order parameter, defined as the 
normalized vector sum of the phase points {Pk(t)) in H{ik{t)— < i >) vs. {ikif)— < i >) 
space, is used to characterize the extent of the synchrony of the population 



where H is the Hilbert transform. This order parameter 
Kuramoto order parameter 
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rZl is similar to the 



10|. The magnitude of r = \Z\, the order, has a maximum value 
of 1 for full synchronization and zero for complete desynchronization (for a population of 
infinite size). 

B. Results 

The population in this study has a nearly unimodal natural frequency distribution (with- 
out coupling and forcing) with a standard deviation of 8.5 mHz; the frequency histogram 
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Figure 1: Experiments: Dimensionless frequency histograms {uj = f / < f > —1) at different 
coupling strengths for a unimodal population. The mean natural frequency < / >= 0.47Hz, 
standard deviation a = 8.5mHz. (a) K = Q. [h) K = 0.026. (c) K = 0.031. (d) K = 0.042. 

is shown in Fig. [1^. Figs. [Tb-d also show the frequency histograms at = 0.026, 0.031 
and 0.042, respectively. Ki K = 0.026 (Fig. [T]b), a dominant, mutually entrained cluster 
emerged at approximately the mean natural frequency. Elements with high and low natural 
frequencies were not entrained yet. With a stronger coupling of = 0.031 (Fig. [lb), the 
cluster grew in size considerably and only a few elements were left out. Ki K = 0.042 
(Fig. did), 62 out of the 64 elements had been in the same cluster with two elements of lower 
natural frequencies desynchronized. External periodic forcing was applied to all these three 
partially synchronized states. 

First consider the forcing of the most synchronized state aX K = 0.042. Without forcing a 
mutually entrained cluster of 62 elements formed at ui =0.450 Hz. The periodic forcing was 
applied at a higher frequency of uj-p = f-p = 0.47Hz. With forcing amplitude of 6 = 3.3mV 
some elements of the original (mutually entrained) cluster moved out and formed a new 
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Figure 2: Experiments: Periodic forcing on a unimodal population with K = 0.042. Vq = 1.09V. 
Mean frequency before forcing < / >= 0.4497 Hz, forcing frequency /f = 0.47Hz, /p— < / >> 2a. 
Top row: Frequency histograms at different forcing amphtudes. Middle row: Frequencies in presence 
of coupling and forcing versus the corresponding natural frequencies. Upper and lower dashed line 
are /p and < / >, respectively. Bottom row: Time series of order r. a. 6 = 3.3mV. h. b = 4.0mV. 
c. b = 4.6mV. d. 6 = 5.3mV. 



cluster at the forcing frequency (forced entrainment) . In Fig. Eh (frequency histogram, top 
panel) there were two large peaks corresponding to the two clusters. The plot of frequencies- 
vs-natural-frequencies (middle panel) also clearly shows two plateaus; the one of the forced 
entrainment was at fp while the other one of the mutual entrainment had a frequency that 
is slightly higher than the mean frequency without forcing. The order oscillated in large 
amplitudes at b = 3.3mV (bottom panel). As b increased to 4.0mV, the forced entrained 
cluster grew in size; simultaneously two mutual entrained clusters emerged (Fig. [2b). The 
frequency of the old mutual entrained cluster (cui) further increased to be closer toup. A new, 
small resonant cluster formed at U2 = loji — up. With a even stronger forcing of 6 = 4.6mV, 
more resonant clusters formed at frequencies of Un = nu\ — {n — l)ujp, n = 2,3,4 (Fig. [2t). 
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Figure 3: Periodic forcing on a unimodal population with weak coupling strengths. Frequencies in 
presence of coupling and forcing versus the corresponding natural frequencies. Top row: K = 0.031, 
< / >= 0.466Hz, /f = 0.486Hz. a. 6 = 3.3mV. h. b = 4.0mV. c. b = 4.6mV. d. 6 = 5.3mV. Bottom 
row: K = 0.026, < / >= 0.4571Hz, /f = 0.48Hz. e. b = 3.3mV. f. 6 = 4.95mV. g. b = 5.3mV. h. 
b = 5.95mV. Upper and lower dashed line are /f and < / >, respectively. 

However, each of these clusters had a smaller size compared with the two mutual entrained 
clusters at b = 4.0mV (Fig. [2)d). With only one large cluster the order no longer exhibited 
large amplitude oscillations. Note that the spacing between any adjacent clusters became 
smaller as b increased. When the forcing amplitude was further increased to b = 5.3mV, 
no distinct resonant clusters were observed and only one large, the forced entrained cluster, 
appeared with 12 non-entrained elements scattering in the lower frequency region (Fig. [2ld). 
Finally, with strong enough forcing amplitude, all the elements formed one cluster at the 
forcing frequency (not shown) . 

During the periodic forcing of a less synchronized state obtained at = 0.031 (Fig. [T]:), 
the occurrence of resonant clusters requires stronger forcing strengths as shown in Figs.[3h-d. 
Forced entrainment occurred with forcing and two main clusters coexisted at weak forcing 
strengths (Figs. [3h, b). As 6 increases, the forced entrained cluster grew larger (Figs.[3h-d) 
and finally became the only dominant cluster at 6 > 5.3mV. The third, new resonant cluster 
with the frequency of uj2 = 1uJ\ — uj^ appeared at 6 = 4.0mV; however, compared with the 
case of = 0.042 at the same weak forcing strength (Fig. [2]P), the new cluster was much 
smaller and can be barely seen (Fig. [3]b). The presence of the third cluster only became 
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obvious at a greater forcing amplitude of 6 = 4.6mV (Fig. [3t); in contrast, in the case of 
K = 0.042 at 6 = 4.6mV multiple resonant clusters had emerged (Fig. [2t). Because of the 
small size of the resonant cluster, the existence of two main clusters at 6 = 4.6mV the order 
oscillated in large amplitudes (not shown). At & = 5.3mV, the forcing was strong enough 
to induce multiple resonant clusters with closer equal frequency spacing (Fig. [3ld) and with 
smaller sizes. Again, before the system reached a fully forced entrainment state (not shown 
in the figures) the mutual entrainment clusters completely disappeared by giving some of 
their elements to the forced entrainment cluster while leaving the others nonentrained. 

The forcing experiments were carried out for an even less synchronized population of 
oscillators obtained with K = 0.026 (Fig. Hb). With a small forcing strength of b = 3.3 
mV (Fig. [3b), a high peak arose at the forcing frequency in the frequency histogram ; 
however, the mutual entrained cluster had broken up because of the weak coupling. In 
the frequency-vs-natural-frequency plot below the forced entrained plateau, the frequencies 
of the nonentrained elements formed an almost continuous 45° line in the lower frequency 
region with a small shoulder around 0.4571 Hz (which is close to the mean natural frequency 
of the oscillators). The order exhibited fairly large oscillations but these oscillations were not 
regular (not shown). As b was increased to a quite strong value of 4.95mV, the elements with 
lower frequencies started to form multiple small groups as shown in Fig. [3]P. However, these 
'clusters' were small and close to each other on the 45° line in the frequency-vs-natural- 
frequency plot (Fig. [3]P). Small clusters were also observed at stronger forcing strengths 
b = 5.3mV (Fig. [3^ ) and b = 5.95mV (Fig. [3h). Since there was no second large cluster 
besides the forced entrained one, the order oscillations were irregular with relatively small 
amplitudes (not shown). Finally, a completely forced entrained cluster was formed with 
strong enough forcing amplitudes. 



III. NUMERICAL SIMULATIONS 



A. Model Equations 

We used a model of anodic electrodissolution of a single nickel electrode proposed by Haim 
et al. [l^. The model in a dimensionless form involves two variables: the dimensionless 
double layer potential drop (e) and the surface coverage of NiO+NiOH (6). One oscillator 
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is described by the following two equations: 

de V — e 



dt R 



1^(9,6) (2) 



^d9 _ exp(0.5e) BChexp{2e) ^ 

dt 1 + Chexp(e) cCh + exp(e) 

where V is the dimensionless applied potential, R is the dimensionless series resistance, T is 
the surface capacity, and ip is the Faradaic current: 



The parameter values = 1600, a = 0.3, B = 6 x 10^^, c = 1 x 10^^ were optimized 
[isl to obtain dynamical features similar to experiments. 

In the exp eriments the oscillators are inherently different because of electrode hetero- 



geneities [l9| and because of addition of different individual external resistors. We model 
the non-identical nature of the oscillators by giving the oscillators different values for the 
parameters R and T in Eqs. [2]and[3l (Other choices are also possible, however, these param- 



eters best approximate the experiments. 
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17l |) For element k the resistance and the surface 
capacitance are obtained using the relationships Rk = {I + ^k)Ro and = (1 + Afc)Fo, 
where is a heterogeneity parameter and Rq and Fq are the mean values. We used a fixed 
value of Fo = 0.01, and Rq = 20 throughout this study. We choose a Lorentzian distribution 
for Afc. (For comparison, simulations with global coupling were also made with a Gaussian 



distribution. |lj|) The Lorentzian distribution p{x) = 7/{7r[(x — XqY + 7^]} is characterized 
by a parameter 7; 27 is the half-width of the distribution. For a typical value of 7 = 0.5, 
the equation parameters Rk and F^ vary within a range of 5% of their means. 

Electrical global coupling of strength K is considered; the model for the coupled set of 



N oscillators is then 
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17|: 



~ir = ~~5~^ " ^F,fc(6'fc, efc) + ^i^(emean(i) " Cfc) (5) 

CLt rti^ iiQ 

d9k exp(0.5efc) 5Chexp(2efc) 
dt l + Chexp(eA;) cGh + exp(efc) 

where Rq is the mean resistance. Global coupling occurs because of the presence of mean 
potential {emean{t) = '^/NYJk=i^k{t)) in Eq. [5l -ft' = represents uncoupled oscillators; 
K ^ CO yields maximum global coupling that synchronizes the oscillators. Global forcing 
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is added through the circuit potential: 



V{t) = Vo + bsm{27ifpt), (7) 

where Vq = 15 is a set potential, b and /f = 0.0745 are the forcing amplitude and frequency, 
respectively. A variable step size fourth order Runge-Kutta method of Matlab was used to 
integrate Eqs. [5] and [6] with a display step size of At = 0.5, absolute error = 10^^, relative 
error = 10^^. Smaller step sizes and error limits gave the same results. A transient of 
t = 10000 was discarded from each time series. 

The phases and frequencies of simulated oscillators are determined from time series data 



of efc(t) with the Hilbert transform method|l5|. 



B. Results 

Without global coupling and forcing the oscillators exhibit frequencies determined by 
the introduced heterogeneities, as shown in Fig. [4^. Since the heterogeneities are small, the 
frequencies of the oscillators follow almost linearly the heterogeneities and thus the frequency 



distribution is similar to a Lorentzian distribution. With global coupling only, a mutual 
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entrained state occurs (see Fig. [Dd) above a critical coupling strength K, 
in this mutually entrained state a large fraction of the oscillators are synchronized. When 
the coupling strength is increased above K^, the number of mutually entrained oscillators 
increases until all the 64 oscillators are fully synchronized. With global forcing only, (see 
Fig. Sb), as the forcing amplitude increases more and more oscillators lock to the forced 
entrainment state. Note however, that there is no phase transition, and the increase in 
the number of entrained oscillators with forcing amplitude is determined by the frequency 
distribution; the large increase at around b = 0.6 is due to the large number of oscillators of 
similar frequencies (peak in Fig. [4^). 

With both global coupling and forcing, in a large fraction of the parameter space {K, b), 
each oscillator can be classified as belonging to one of the three major groups: desynchronized 
(D), mutually entrained (M), and forced entrained (F) states. Fig. [Hd shows the regions in 
the {K, b) parameter space in which these states can occur. At low forcing and coupling 
the system is desynchronized (D). Starting from this desynchronized state with increasing 
coupling there is a transition to the occurrence of the mutually entrained state (M+D in 
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Figure 4: (Color Online) Numerical Simulations: Mutual and forced entrainment in a model of sixty- 
four globally coupled oscillators with external forcing, a. Natural (inherent) frequency distribution, 
b. Mutual entrainment with global coupling only. Number of mutually entrained oscillators vs. 
coupling strength K. c. Forced entrainment with external forcing only. Number of oscillators 
locked to external forcing signal vs. forcing amplitude b. d. Phase diagram showing forced and 
mutually entrained regions in the coupling strength — forcing amplitude parameter space. F: Forced 
entrainment. M: Mutual entrainment. D: desynchronized state, fp = 0.0745. 

Fig. [Ub). With increasing the forcing amplitude in a weakly coupled population, a forced 
entrainment state starts to occur (F+D). At very strong coupling, the transition is similar 
to the forcing of a single oscillator: the frequency of the mutually entrained cluster (M) 
shifts to that of the forcing (F). 

We investigated the dynamics in the middle region of the phase diagram further, where 
both forced and mutually entrained states can co-exist. (Note that this is not hysteresis. In 
a synchronized state some of the oscillators are mutually, some are forced entrained; there 



10 



50 




1. 



















1 . 






72 


74 76 78 
10^ f 


72 74 76 78 
10^ f 


72 74 76 78 
10^ f 


78 
'^-^ 76 

CO 

o 74 




— — — 1-. — i- — - 


m 

m * 

• • 


— • — *- -•- -«-««M 


• 

• 

• 

-:' 


72 













70 „ 72 74 
10^ r(ib = 0) 



70 „72 74 70 „ 72 74 

^o^f.{b = o) f.{b = o) 











10"^ t 




10 



10"^ t 



10 



10 



10"^ t 



Figure 5: Numerical Simulations: frequency distributions and order parameters at K = 0.67. 
Top row: frequency distribution of oscillators. Middle row: frequency of oscillators vs. inherent 
frequency (dashed line denotes the forcedly entrained cluster). Bottom row: order parameter vs. 
time. a. Mutual entrainment dominated state with weak forcing, b = 0.428. b. Resonance 
clustering at intermediate forcing strength, b = 0.443. c. Entrainment to strong forcing, b = 0.452. 



are some oscillators that are also entrained but, as we shall see below, their frequencies are 
combination of those of the mutualy and forced entrainment states.) The frequencies and 
the order of such a state with weak forcing is shown in Fig. [5h. Since the forcing is weak, 
the large peak in the frequency histogram is associated with the mutually entrained state, 
while the smaller peak at the forcing frequency with the forced entrainment. The frequency 
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vs. inherent frequency graph (middle) shows that the elements with lower frequencies form 
the mutually entrained state; such an asymmetry in the entrainment was also observed with 
pure coupling as well, [l^ Elements with frequencies between the mutually and forced 
entrainment groups have strongly scattered, lower frequencies; those oscillators with larger 
frequencies are not entrained. Because the majority of the population is mutually entrained, 
a large value of order parameter (bottom panel) is observed with small modulation due to 
forced entrainment state. 

At a somewhat stronger forcing amplitude, a restructuring of the frequencies of the oscil- 
lators occurs (see Fig. [5]b). Now the major group is associated with the forced entrainment, 
however, the mutually entrained state is difficult to determine since the population becomes 
clustered at frequencies with a (dimensionless) spacing of approximately 0.0009 [/ ={0.0736 
, 0.0745 (forcing), 0.0754 , 0.0764]; 8 elements of the population of 64 oscillators do not be- 
long to these resonant clusters. Resonant clusters seem to appear as an interaction between 
the mutual and forced entrainments and produce a strongly oscillating order parameter. 

As the forcing amplitude is further increased, the forced entrainment (Fig. Eb) state 
dominates in which the elements around the forcing frequencies along with those of the lower 
frequencies are entrained and the elements with higher frequencies are not synchronized. The 
order parameter again has a large, close to 1 value, but with small fluctuations due to the 
desynchronized elements. 

Resonance clustering to a lesser extent also appears at weaker coupling strength as it is 
shown at three forcing amplitudes in Fig. [6l At these conditions (Figs. Ehb), the elements 
do form various frequency clusters, however, there may be quite a few clusters with smaller 
number of elements and the spacing between the clusters is small. These states also exhibit 
oscillating order (not shown in the figures). At a state very close the destruction of mutual 
entrainment (Fig. [Bh) the frequencies of mutual and forced entrainment states are very 
close and frequency clustering appears as discretized frequencies of groups of 1-3 oscillators 
only. 

Fig. [7] shows the simulation results with a population of 512 oscillators at inherent 
frequency distribution, coupling strength, and forcing amplitude similar to those investigated 
with 64 oscillators. We observed five frequency clusters with approximate spacing of 0.001 
[at frequencies /={0.0734, 0.0745 (forcing), 0.0756, 0.0766, 0.0777}]. This 5-cluster state 
with = 512 is reminiscent of the 4-cluster state with A^ = 64 (Fig. EJd), however, one 
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Figure 6: Numerical simulations: frequencies of oscillators (with coupling and forcing) vs. inherent 
frequencies at weaker interactions, K = 0.54. a. Weak forcing, b = 0.385. b. Intermediate forcing 
strength, b = 0.398. c. Strong forcing, b = 0.410. 
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Figure 7: Numerical simulations: population of 512 oscillators exhibiting resonance clustering at 
K = 0.67 and b =0.47. a. Frequency distribution of oscillators. (For / = 0.075 there is a large 
peak with N = 325, therefore, a zoom is shown) b. Frequency of oscillators vs. inherent frequency 
(dashed line denotes the forcing frequency). 

additional cluster could be resolved at / = 0.0777. Thus, the fine structure of resonance 
clusters can be better seen with larger populations. 



IV. DISCUSSION 



Clustering is a behavior of a population of oscillators through which dynamically differen- 
tiated elements form. The classification of clusters [2]Jis based on whether it affects mainly 
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the phases (phase clusters) [22, l23| or the amplitudes (amplitude clusters) [2j| of the peri- 



odic oscillators. Chaotic systems are particularly inclined for producirig c 



usters at coupling 

strengths weaker than those required for identical synchronization[25|, In this study, 

we observed frequency clusters in a forced, weakly coupled population of limit-cycle electro- 
chemical oscillators with unimodal heterogeneities. In addition to the trivial mutually and 
entrained clusters, groups of elements are obtained at other, discrete frequency values. We 
also showed that at coupling strengths and forcing amplitudes with well-defined resonance 
clusters the frequencies are equally spaced and follow the relationship found by Sakaguchi 
lll|: Un = nuJi — (n — 1)u!f. (Such spectrum was analyzed by linear stability of the forced 



entrainment state. 



28l |) With weaker coupling strengths a large number of clusters with small 



number of elements were observed whose frequency was not greatly modified from the nat- 
ural frequencies. We note that frequency discretization through a different mechanism, viz., 



through large delays in the coupling term, was also observed in a pair of lasers 29jand in a 
pair chaotic Rossler oscillators 30|. 

Along with the occurrence of frequency clusters, a strongly oscillating order parameter 
was observed; the order decreased to very low values and thus intermittent loss of the overall 
rhythm could be observed. The loss of the overall rhythm may have implications in the 
dynamics of biological systems where the rhythm could be either essential or pathological 

i 

A biological example in which coupling and forcing play a role is the circadian master 
clock in the brain[3]J; the suprachiasmatic nuclei consist of heterogeneous, circadian oscil- 
lators that are entrained by cell interactions and by external light. The SCN is strongly 
heterogeneous and the interaction of the core, being entrained by light, and the shell (com- 
posed of mostly self-oscillating neurons) is a complex issue|32|,|33|; resonance clustering could 
play a role in the dynamical behavior of the different regions. Resonant interactions between 



oscillators 
tory units 



lave been considered as well in the context of information processing with oscilla- 
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the cortex is considered as a weakly coupled oscillator population forced 
by the thalamic input. 

As far as the importance of resonance clustering in biological networks of rhythmic el- 
ements is concerned, we shall consider that we investigated a global coupling topology. 
Although such global, all-to-all coupling is not very likely in biological systems, sometimes 
a network of oscillators can be approximated by global coupling. A fundamental communi- 
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cation mechanism of bacteria, quorum sensing, is often modeled by global interactions; such 
a mechanism was shown to be able to produce synchrony js^. Synchronization of circadian 
cells was modeled by a global interaction mechanism based on an argument that the spatial 
transmission of neurotransmitters released by each cell is fast compared to the timescale 



of oscillations [37|. Neural networks with electrically spiking neurons are also often can be 



considered as a population of weakly, coupled globally coupled oscillators 



Even when 



network coupling occurs, the depth of the network is often not very large when external 



forcing is effective [39| thus global coupling approximation and resonance clustering can play 
a role in generation of collective dynamics with coupling and forcing. 

Some features of the cluster interactions such as generation of the resonant clusters could 
also be observed in a population with bimodal heterogeneities; the forced system can be re- 
garded as a simplified case of bimodal system where one special mode has standard deviation 
of zero and its strength can be directly controlled. Frequency clustering and generation of 
complex collective signal were found with the global coupling of periodic electrochemical os- 
cillator populations with bimodal natural frequency distributions 4^; merging and splitting 
of clusters occurred on the way to the final synchronized state with increasing the coupling 
strength. Such bimodal (and multimodal) populations may occur in biological systems that 
are composed of broadly heterogeneous cell groups. 

Although we investigated resonance clustering here with electrochemical oscillators, sim- 
ilar dynamical differentiation mechanisms are expected to occur in a variety of rhythmic 
multicellular systems under the co-operative effects of coupling and forcing; the synchro- 
nized groups of elements would contribute to the formation of multi-structured hierarchical 



organizations often seen in complex systems 



Acknowledgments 



41 



421 



Financial support from the National Science Foundation (CBET-0730597) is acknowl- 
edged. 



[1] A. T. Winfree, The geometry of biological time (Springer- Verlag, New York, 1980). 
[2] V. Petrov, Q. Ouyang, and H. L. Swinney, Nature 388, 655 (1997). 



15 



[3] A. von Oertzen, H. H. Rotermund, A. S. Mikhailov, and G. Ertl, J. Phys. Chem. B 104, 3155 
(2000). 

[4] A. L. Lin, M. Bertram, K. Martinez, H. L. Swinney, A. Ardelea, and G. F. Carey, Phys. Rev. 
Lett. 84, 4240 (2000). 

[5] B. Marts, A. Hagberg, E. Meron, and A. L. Lin, Phys. Rev. Lett. 93, 108305 (2004). 

[6] V. K. Vanag, A. M. Zhabotinsky, and L R. Epstein, Phys. Rev. Lett. 86, 552 (2001). 

[7] O. Steinbock, V. Zykov, and S. C. Miiller, Nature 366, 322 (1993). 

[8] O. Rudzick and A. S. Mikhailov, Phys. Rev. Lett. 96, 018302 (2006). 

[9] H. Fukuda, N. Tamari, H. Morimura, and S. Kai, J. Phys. Chem. A 109, 11250 (2005). 
[10] Y. Kuramoto, Chemical oscillations, waves and turbulence (Springer, Berhn, 1984). 
[11] H. Sakaguchi, Prog. Theoret. Phys. 79, 39 (1988). 

[12] E. Montbrio, J. Kurths, and B. Blasius, Phys. Rev. E 70, 056125 (2004). 

[13] L Z. Kiss, Y. M. Zhai, and J. L. Hudson, Science 296, 1676 (2002). 

[14] L Z. Kiss, W. Wang, and J. L. Hudson, J. Phys. Chem. B 103, 11433 (1999). 

[15] A. S. Pikovsky, M. G. Rosenblum, and J. Kurths, Synchronization: A Universal Concept in 

Nonlinear Science (University Press, Cambridge, 2001). 
[16] Y. M. Zhai, L Z. Kiss, and J. L. Hudson, Ind. Eng. Chem. Res. 43, 315 (2004). 
[17] Y. M. Zhai, L Z. Kiss, H. Daido, and J. L. Hudson, Physica D 199, 387 (2004). 
[18] D. Haim, O. Lev, L. M. Pismen, and M. Sheintuch, J. Phys. Chem. 96, 2676 (1992). 
[19] L Z. Kiss, W. Wang, and J. L. Hudson, J. Phys. Chem. B 103, 11433 (1999). 
[20] A. T. Winfree, J. Theor. Bioh 16, 15 (1967). 
[21] A. S. Mikhailov and K. Showalter, Phys. Rep. 425, 79 (2006). 

[22] D. Golomb, D. Hansel, B. Shraiman, and H. Sompolinsky, Phys. Rev. A 45, 3516 (1992). 
[23] K. Okuda, Physica D 63, 424 (1993). 

[24] V. Hakim and W. J. Rappel, Phys. Rev. A 46, R7347 (1992). 
[25] K. Kaneko, Phys. Rev. Lett. 63, 219 (1989). 

[26] K. Kaneko and L Tsuda, Complex Systems: Chaos and Beyond (Springer, Berlin, 2001). 
[27] W. Wang, L Z. Kiss, and J. L. Hudson, Chaos 10, 248 (2000). 

[28] T. M. Antonsen Jr, R. T. Faghih, M. Girvan, E. Ott, and J. Platig, submitted to Chaos 

(|arXiv:0711.4135ll (2008). 
[29] H. J. Wiinsche, S. Bauer, J. Kreissl, O. Ushakov, N. Korneyev, F. Henneberger, E. Wille, 

16 



H. Erzgraber, M. Peil, W. Elsasser, et al., Phys. Rev. Lett. 94, 163901 (2005). 
[30] S. Yanchuk, Phys. Rev. E 72, 036205 (2005). 
[31] S. M. Reppert and D. R. Weaver, Nature 418, 935 (2002). 
[32] M. C. Antle and R. Silver, Trends in Neurosciences 28, 145 (2005). 

[33] W. Nakamura, S. Yamazaki, N. N. Takasu, K. Mishima, and G. D. Block, J. Neurosci. 25, 
5481 (2005). 

[34] F. C. Hoppensteadt and E. M. Izhikevich, Biosystems 48, 85 (1998). 
[35] F. C. Hoppensteadt and E. M. Izhikevich, Phys. Rev. Lett. 82, 2983 (1999). 
[36] J. Garcia-Ojalvo, M. B. Elowitz, and S. H. Strogatz, Proc. Natl. Acad. Sci. USA 101, 10955 
(2004). 

[37] D. Gonze, S. Bernard, C. Waltermann, A. Kramer, and H. Herzel, Biophys. J. 89, 120 (2005). 
[38] H. Kori, Phys. Rev. E 68, 021919 (2003). 

[39] H. Kori and A. S. Mikhailov, Phys. Rev. Lett. 93, 254101 (2004). 

[40] A. S. Mikhailov, D. H. Zanette, Y. M. Zhai, I. Z. Kiss, and J. L. Hudson, Proc. Nat. Acad. 
Sci. USA 101, 10890 (2004). 

[41] S. C. Manrubia, A. S. Mikhailov, and D. H. Zanette, Emergence of Dynamical Order: Syn- 
chronization Phenomena in Complex Systems (World Scientific, Singapore, 2004). 

[42] K. Kaneko, Life: An Introduction to Complex System Biology (Springer, 2006). 



17 



